Viscous fingering of miscible slices 
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Viscous fingering of a miscible high viscosity slice of fluid displaced by a lower viscosity fluid is 
studied in porous media by direct numerical simulations of Darcy's law coupled to the evolution 
equation for the concentration of a solute controlling the viscosity of miscible solutions. In contrast 
with fingering between two semi-infinite regions, fingering of finite slices is a transient phenomenon 
due to the decrease in time of the viscosity ratio across the interface induced by fingering and 
dispersion processes. We show that fingering contributes transiently to the broadening of the peak 
in time by increasing its variance. A quantitative analysis of the asymptotic contribution of fingering 
to this variance is conducted as a function of the four relevant parameters of the problem i.e. the 
log-mobility ratio R, the length of the slice I, the Peclet number Pe and the ratio between transverse 
and axial dispersion coefficients e. Relevance of the results is discussed in relation with transport 
of viscous samples in chromatographic columns and propagation of contaminants in porous media. 



I. INTRODUCTION 

Viscous fingering is an ubiquitous hydrodynamic in- 
stability that occurs as soon as a fluid of given viscosity 
displaces another more viscous one in a porous medium 
0. As such, the typical example usually presented for 
this instability is that of oil recovery for which viscous 
fingering takes place when an aqueous solution displaces 
more viscous oil in underground reservoirs. This explains 
why numerous articles devoted to the theoretical and 
experimental analysis of fingering phenomena have ap- 
peared in the petroleum engineering community [lj. For 
what concerns the geometry, theoretical works typically 
focus on analyzing the stability properties and nonlin- 
ear dynamics of an interface between two semi-infinite 
domains of different viscosity. In the same spirit, exper- 
imental works done either in real porous media or in a 
model Hele-Shaw system (two parallel plates separated 
by a thin gap width) consist in injecting continuously a 
low viscous fluid into the medium initially filled with the 
more viscous one. The attention is then focused on the 
dynamics of the interface between the two regions. The 
instability develops and the fingers grow continuously in 
time until the displacing fluid has invaded the whole ex- 
perimental system. As long as the experiment runs {i.e. 
until the displacing fluid reaches the outlet), the insta- 
bility develops. Dispersion of one fluid into the other 
may lead to a slight stabilization in time nevertheless 
this stabilization is usually negligible on the time scale of 
the experiment and for high injection rates. 

The situation is drastically different in other impor- 
tant applications in which viscous fingering is observed, 
such as in liquid chromatography or groundwater con- 
tamination. Liquid chromatography is used to separate 
the chemical components of a given sample by passing it 
through a porous medium. In some cases, and typically 
in preparative or size exclusion chromatography, the vis- 
cosity of the sample is significantly different than that 
of the displacing fluid (the eluent). Displacement of the 



sample by the eluent of different viscosity leads then to 
viscous fingering of either the front or the rear interface 
of the sample slice, leading to deformation of the initial 
planar interface. This fingering is dramatic for the per- 
formance of the separation technique as it contributes to 
peak broadening and distortions. Such conclusions have 
been drawn by several authors that have shown either 

experimentally HHSSH00 or numericall y 
the influence of viscous fingering on peak deformations. 

In groundwater contamination and more generally soil 
contamination, it is not rare that the spill's extent is finite 
due to a contamination localized in space and/or time. If 
the spill's fluid properties are different than that of water, 
and in particular, if they have different viscosity and/or 
density H G3, fingering phenomena may influence the 
spreading characteristics of the contaminated zone. For 
ecological reasons, it is important then to quantify to 
what extent fingering will enlarge the broadening in time 
of this polluted area. 

Nonlinear simulations of fingering of finite samples 
have been performed in the past by Tucker Norton, Fer- 
nandez et al. 0, @ in the context of chromatographic 
applications, by Christie et al. in relation to "Water- 
Alternate Gas" (WAG) oil recovery techniques ^1 as 
well as by Zimmerman [lj and have shown the influence 
of fingering on the deformation of the sample without 
however investigating the asymptotic dynamics. Man- 
ickam and Homsy, in their theoretical analysis of the sta- 
bility and nonlinear dynamics of viscous fingering of mis- 
cible displacements with nonmonotonic viscosity profiles 
have further stressed the importance of reverse fingering 
in the deformation of finite extent samples [l3L[T^ |. Their 
parametric study has focused on analyzing the influence 
of the endpoint and maximum viscosities on the growth 
rate of the mixing zone. 

In this framework, the objectives of this article are 
twofold: first, we analyze the nonlinear dynamics of vis- 
cous fingering of miscible slices in typical analytical chro- 
matographic and groundwater contamination conditions 
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in order to underline its specificities and, second, we 
quantify the asymptotic contribution of viscous finger- 
ing to the broadening of the output peaks as a function 
of the important parameters of the problem. From a nu- 
merical point of view, the only difference with regard to 
most of the previous works devoted to viscous fingering 
[TEl fill IT?! is the initial condition which is now a sample 
of finite extent instead of the traditional interface be- 
tween two semi-infinite domains. As we show, this has 
an important consequence: if the longitudinal extent of 
the slice is small enough with regard to the length of the 
migration zone, dispersion becomes of crucial importance 
as it leads to such a dilution of the displaced sample into 
the bulk fluid before reaching the measurement location 
that fingering just dies out. As a consequence fingering is 
then only a transient phenomenon and the output peak of 
the diluted sample may look Gaussian even if its variance 
is larger than that of a pure diffusive dynamics because 
of transient fingering. This explains why the importance 
of fingering phenomena in chromatography and soil con- 
tamination has been largely underestimated or ignored in 
the literature. We perform here numerical simulations to 
compute the various moments of the sample distribution 
as a function of time when fingering takes place. This 
allows us to extract the contribution of viscous finger- 
ing to the variance of the averaged concentration profile 
and to understand how this contribution varies with the 
important parameters of the problem which are the log- 
mobility ratio R between the viscosity of the sample and 
that of the bulk fluid, the Peclet number Pe, the dimen- 
sionless longitudinal extent I of the slice and the ratio 
e between the transverse and longitudinal dispersion co- 
efficients. The outline of the article is the following. In 
Sec.m we introduce the model equations of the problem. 
Typical experimental parameters for liquid chromatogra- 
phy and groundwater contamination applications are dis- 
cussed in Sec. IIIII The characteristics of the fingering of 
a miscible slice are outlined in Sec. lIVI while a discussion 
on the moments of transversely averaged profiles is done 
in Sec. [3 Eventually, a parametric study is conducted 
in Sec. IvH before a discussion is made. 



II. MODEL SYSTEM 

Our model system is a two-dimensional porous medium 
of length L x and width L y (Fig. . A slice of fluid 2 of 
length W is injected in the porous medium initially filled 
with carrier fluid 1. This fluid 2, which is a solution of 
a given solute of concentration c-i in the carrier, will be 
referred in the following as the sample. This sample is 
displaced by the carrier fluid 1 in which the solute concen- 
tration c is equal to c\ = 0. Assuming that the viscosity 
of the medium is a function of the concentration c and 
that the flow is governed by Darcy's law, the evolution 
equations for the system are then: 




FIG. 1: Sketch of the system. 
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where /i is the viscosity of the fluid, K is the permeability 
of the medium, p is the pressure and u — (u, v) is the two- 
dimensional velocity field. The displacing fluid is injected 
in a uniform manner with a mean velocity U along the 
x direction. D x ,D y are the dispersion coefficients along 
the flow direction and perpendicular to it respectively. 
The characteristic speed U is used to define a charac- 
teristic length L c = D x /U and time r c = D x /U 2 . We 
nondimcnsionalize space, speed and time by L c , U and 
r c respectively. Pressure, viscosity and concentration are 
scaled by niD x /K 7 /xi and C2, where fix is the viscosity 
of the displacing fluid and C2 the initial concentration of 
the sample. The dimensionless equations of the system 
become 



dc 
di 



V-M = 

Vp = 
u ■ Vc — 



o, 

-H(c)u, 
d 2 c 
dx 2 



d^c 
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where e = D y /D x . If e = 1, dispersion is isotropic while 
e 7^ 1 characterizes anisotropic dispersion. Switching to 
a coordinate system moving with speed U, i.e. making 
the change of variables x' — x — t, y = y, u' = u — i x with 
i x being the unit vector along x, we get, after dropping 
the primes: 



dc 
di 



Vp = 
u • Vc = 



0, 



-/i(c)(u + i x ), 

d 2 c d 2 c 
dx 2 dy 2 ' 



(7) 

(8) 

(9) 



We suppose here that the viscosity is an exponential func- 
tion of c such as: 



/x(c) 



(10) 



0. 



(1) 



where R is the log-mobility ratio defined by R = 
Zn(/i2 //•*].), where \ii is the viscosity of the sample and, 
as said before, /ii is the viscosity of the displacing fluid 
(Fig. 0. If R > 0, then we have a low viscosity fluid dis- 
placing a high viscosity sample and the rear interface of 
the sample will be unstable with regard to viscous finger- 
ing. If R < 0, then the sample is the less viscous fluid and 
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the front interface of the slice will then develop fingering. 
In our simulations, we consider the R > situation. 

Introducing the stream function ip such that u = 
dip/dy and v — —dip/dx^ taking the curl of Eq. JHJl, we 
get our final equations |18| : 



V 2 ip = R 

dc 
at 



dip dc dip dc 
dx dx dy dy 
dtp dc dip dc d 2 c 
dy dx dx dy dx 2 



dy 



_c?c 
dy 2 



(11) 
(12) 



This model is numerically integrated using a pseudo- 
spectral code introduced by Tan and Homsy and 
successfully imp lemented for various numerical studies 
of fingering |19l |20| . The two-dimensional domain of in- 
tegration is, in dimensionless units, of size Pe x L where 
Pe = ULy/ D x is the dimensionless width which is noth- 
ing else than the Peclet number of the problem, while 
L = UL X /D X . The dimensionless length of the sam- 
ple is I = UW/D X . The initial condition corresponds to 
a convectionless fluid ("0 = everywhere) embedding a 
rectangular sample of concentration c = 1 and of size 
Pe x I in a c = background. The middle of the sample 
is initially located at x = 2L/3. In practice, for the sim- 
ulations, the initial condition corresponds to two back to 
back step functions between c = and c = 1 with an 
intermediate point where c — h+A-r, r being a random 
number between and 1 and A the amplitude of the noise 
(typically of the order of 10~ 3 ). This noise is necessary 
to trigger the fingering instability on reasonable comput- 
ing time. If A = 0, numerical noise will ultimately seed 
the fingering instability but on a much longer time scale. 
The boundary conditions are periodic in both directions. 
This is quite standard along the transversal direction y. 
This does not make any problem along the x-axis as c = 
at both x = and x = L. The problem is controlled by 
four dimensionless parameters: the log-mobility ratio R, 
the Peclet number Pe, the initial length of the injected 
sample I and the ratio between transverse and longitudi- 
nal dispersion coefficients e. 



III. EXPERIMENTAL VALUES OF 
PARAMETERS FOR TWO APPLICATIONS 

In order to perform numerical simulations, let us com- 
pute the order of magnitude of the main parameters 
(Peclet number Pe, length of the sample I) for both a 
liquid chromatography experiment and for the propaga- 
tion of contaminants in a porous medium (groundwater 
contamination) . 



A. Chromatographic applications 

First of all, let us note that in most chromatographic 
applications, heterogeneous chemistry (particularly ad- 
sorption and desorption phenomena) is crucial to the 



separation process and will undoubtedly affect possible 
fingering processes. We neglect such physicochemical in- 
teractions in this first approach focusing on the effect of 
viscous fingering on an unretained compound. A typical 
chromatographic column has a diameter d = 4.6 mm, a 
length L x — 150 mm and consists of a porous medium 
packed with porous particles, the total (intraparticle and 
interparticle) porosity being equal to 0.7. The volume of 
the sample introduced in the column is of the order of 
20/ul, injected at a flow rate Q ~ 1ml min -1 . The ex- 
tent of the injected sample is then W ~ 1.7 10 -3 m and 
the speed of the flow U ~ 1.4 10~ 3 ms -1 . The longitu- 
dinal and transverse dispersion coefficients are typically 
D x = 1.43 lO^n^s" 1 (22 and D y = 5.65 10~ 10 m 2 s" 1 
(22| . These parameters allow one to define a character- 
istic length L c — D x /U and a characteristic time r c = 
D x /U 2 . As a result, the Peclet number Pe — Ud/D x 
is here nothing else than the dimensionless diameter i.e. 
Pe = d/L c ~ 460, while the dimensionless longitudinal 
extent of the sample becomes I = W/L c ~ 170. The dis- 
persion ratio is equal to e = D y /D x ~ 0.04. As a typical 
transit time from inlet to outlet takes roughly r = 100 s, 
the dimensionless time of a simulation should be of the 
order of T = t/t c ~ 15000 to account for a realistic time 
to characterize the properties of the output peaks. 



Soil contamination 



The effects of fluid viscosity and fluid density may 
be important in controlling groundwater flow and so- 
lute transport processes. Recently, a series of column 
experiments were conducted and analyzed by Wood et 
al. 01 to provide some insight into these questions. The 
experiments were performed in fully saturated, homo- 
geneous and isotropic sand columns (porosity equals to 
0.34 and e = 1) by injecting a 250 ml pulse of a known 
concentration solution at a flow rate Q — 0.015 m 3 /day. 
Their experimental setup consists of a vertical pipe L x = 
0.91m in length with a diameter of d = 0.15 m. As- 
suming the medium to be homogeneous and the disper- 
sion coefficient D as isotropic, a typical value for the 
aquifer dispersion coefficient is D = 0.1m 2 /month i.e. 
D ~ 3.86 10~ 8 m 2 s _1 (2^. In the same spirit as above, 
we compute the Peclet number to be of the order of 
Pe ~ 110, while the dimensionless length of the sample 
is / ~ 30. 

Based on these two examples, let us now investigate the 
properties of fingering of finite slices for typical values 
of parameters in the range computed above i.e. Pe ~ 
100 - 500, e ~ 0.04 - 1, I ~ - 500, while R is supposed 
to be of order one. 



FIG. 2: Density plots of concentration at successive times in 
the frame moving at the velocity U. From top to bottom: t~0, 
500, 700, 1000, 1500, 2000, 5000, 15000 and 60000 (Pe = 512, 
1= 128, ii = 2,e = l). 



IV. FINGERING OF A FINITE WIDTH 
SAMPLE 

Figure [5] shows in a frame moving with the injection 
velocity U the typical viscous fingering of the rear inter- 
face of a sample displaced from left to right by a less vis- 
cous fluid. The system is shown at successive times using 
density plots of concentration with black (resp. white) 
corresponding to c = 1 (resp. c = 0). While the front in- 
terface is stable, the back interface develops fingers such 
that the center of gravity of the sample is displaced in the 
course of time towards the back with regard to its initial 
position. This dynamics results from the fact that the 
stable zone acts as a barrier to finger propagation in the 
flow direction leading therefore to reverse fingering. Such 
a reverse fingering has been well characterized by Man- 
ickam and Homsy in their numerical analysis of fingering 
of nonmonotonic viscosity profiles |l4j . After a while, 
dispersion comes into play and dilutes the more viscous 
fluid into the bulk of the displacing fluid. As the sample 
becomes more and more diluted, the effective viscosity 
ratio decreases in time weakening the source of the in- 
stability. Ultimately, dispersion becomes dominant and 
the sample goes on diluting in the bulk without witness- 
ing any further fingering phenomenon. These successive 
steps can clearly be observed on the transverse averaged 



FIG. 3: Transverse average profiles of concentration at suc- 
cessive times t=0, 500, 700, 1000, 1500, 2000, 5000 and 
15000. Insert: transverse average profile of concentration at 
t = 60000 (Pe = 512, I = 128, R = 2, e = 1). 



profiles of concentration defined as 
1 f Pe 

c ( a: '*) = p^y c(x,y,t)dy. (13) 

As seen on Fig. we first start with two back to back 
step functions defining an initial sample of extent I. Dur- 
ing the first stages of the injection, there is a first diffu- 
sive regime quickly followed by the fingering of the rear 
interface corresponding here to the left front. While the 
right (i.e. frontal) interface features the standard error 
function characteristic of simple dispersion, the left one 
shows bumps signaling the presence of fingering. Because 
the extent of the sample is finite, dispersion and fingering 
contribute to the fact that the maximum concentration 
becomes smaller than one, effectively leading to a viscos- 
ity ratio between the sample and the bulk that decreases 
in time. As a consequence, fingering dies out and the 
transverse profile starts to follow a distorted Gaussian 
shape. If one waits long enough, the asymmetry of the 
bell shape diminishes which explains that output peaks 
in chromatographic columns may look Gaussian even if 
fingering has occurred during the first stages of the travel 
of the sample in the column. As computed in the pre- 
ceding section, a typical dimensionless time of transit 
in a real chromatographic setup corresponds to 15000 
units of time. Figures |3 and show that, after 15000 
units of time, fingering is disappearing for this specific 
set of typical values of parameters, and that dispersion 
becomes again the dominant mode. As chromatographic 
columns are generally opaque porous media, it is there- 
fore not astonishing that the presence of viscous fingering 
has long been totally ignored until recent experimental 
works which have visualized fingeringby magnetic reso- 
nance or optical imaging US USES 11 Similarly, 
tracing of the spatial extent of a contaminant plume at a 
distance far from the pollution site may lead to measure- 
ments of Gaussian-type spreading even if fingering has 
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occurred at early times. The only influence of such fin- 
gering appears in the larger variance of the sample than 
in the case of pure dispersion as we show it next. 



V. MOMENTS OF THE TRANSVERSE 
AVERAGED PROFILE 

The averaged profiles of concentration c(x, t) allow to 
compute the three first moments of the distribution: the 
first moment m 



M Jo c ( x > l ) xdx 
m(t) = — — 

J c(x, t)dx 



(14) 



is the position of the center of mass of the distribution as 
a function of time. The second moment is the variance 



J Q c(x,t)[x — m(t)] 2 dx 
c(x, t)dx 



a 2 (t) 



(15) 



giving information on the width of the distribution. 
Eventually, we compute also the third moment, i.e. the 
skewness 



a(t) = 



c(x,t)[x — m(t)] 3 dx 
c(x, t)dx 



(16) 



that gives information concerning the asymmetry of the 
peak with regard to its mean position. 

The variance a 2 is the sum of three contributions: 



a 2 {t) 
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(17) 



where a 2 = I 2 /12 is the variance due to the initial length 
of the sample, a 2 , = 2t is the contribution of dispersion in 
dimensionless units and a 2 is the contribution due to the 
fingering phenomenon. If R = 0, the displacing fluid and 
the sample have the same viscosity and no fingering takes 
place. Hence, in that case, a 2 = o 2 + a 2 d = l 2 /l2 + 2t. We 
have checked that this result is recovered by numerical 
simulations for R = 0. The integrals in the computation 
of the moments I|14[) -(I16 [1 are evaluated numerically by 
using Simpson's rule. The numerical result is very good if 
the spatial discretization step dx is small. Typically, we 
get the exact result for dx = 1. Unfortunately, dx = 1 is 
a resolution too high for fingering simulations especially 
if one wants to look at the dynamics at very long times. 
As an example, pre vious simulations on viscous fingering 
phenomena |15l Il9j were done with larger dx as typical 
dimensionless fingering wavelengths are around 100 for 
R = 3 for instance. Using typically dx = 4 gives roughly 
25 points per wavelength which is numerically reasonable. 
For what concerns the variance, simulations with R = 
and dx = 4 give the correct of at t — but a constant 
shift appears so that cr 2 (t) — I 2 / 12 — 2t = C, with C 
being a constant of the order of 0.1% of l 2 /12. As we are 




-60x 1 (P 



50 60x1 1 



FIG. 4: First three moments of the distribution: (a) Mean 
position m of the center of mass, (b) Variance a 2 , (c) Skewness 
a (Pe = 512, I = 128, R= 2, e = 1). 



mostly interested in the rate of variation of cr/, where a/ 
is defined as 



«f{t) = Jo}(t) = Ja 2 (t) 



(18) 



all simulations are done here with dx — 4. The slight C 
shift does not affect the value of cr/ as we have checked 
it for decreasing values of dx. 

Figurc0|shows the temporal evolution of the first three 
moments i.e. the deviation m(t) — m(t = 0) of the center 
of mass in comparison to its initial location at t — 0, the 
total variance <x 2 (i) and the skewness a(t) for the typical 
example of Figs. |3 and As fingering occurs quicker 
than dispersion, the center of gravity of the sample m(t) 
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r 



FIG. 5: Mixing zone Ld as a function of time realized with 
the same parameters (Pe = 512, I = 128, R = 2, e = 1) 
and three different values of the amplitude A of the noise 
seeding the initial condition: (• ■ •) A = 0.1, (- -) A = 0.01, 
( — ) A — 0.001; (- • -) theoretical curve of a pure diffusive 
behavior Ld oc \ft. Insert: zoom on the first stages on the 
injection. The onset time t* , corresponding to the time at 
which the mixing zone departs from the pure diffusive initial 
transient, is a decreasing function of A. 



is displaced towards the back (smaller x values) because 
of reverse fingering of the rear interface of the sample 
[Fig. @Ja)]. Fingering contributes to the widening of the 
peak and thus a 2 (t) increases [Fig.^Jb)] while the skew- 
ness a(t) becomes non zero due to the asymmetry of the 
fingering instability with regard to the middle of the sam- 
ple [Fig.0Jc)]. After a while, fingering dies out and the 
first moment m(t) saturates to a constant indicating that 
dispersion becomes again the only important dynamical 
transport mechanism. Note that the skewness a is ob- 
served to revert back towards at very long times. 

Onset of fingering is also witnessed in the growth of 
the mixing zone Ld defined here as the interval in which 
c(x,t) > 0.01 (Fig. ISJ. An important thing to note is 
that, after a diffusive transient, fingering appears on a 
characteristic time scale t* , defined as the time for which 
the mixing zone temporal dependence departs from the 
pure diffusive regime. 

As has already been discussed before the char- 

acteristics of the fingering onset time t* and of the de- 
tails of the nonlinear fingering regime are dependent on 
the noise amplitude A. The higher the noise intensity A, 
the quicker the onset of the instability (Insert in Fig.JSJ. 
To get insight into the influence of the relevant physical 
parameters of the problem, it is therefore necessary to 
fix the amplitude of the noise to an arbitrary constant 
as this is not a variable that is straightforwardly experi- 
mentally available. In that respect, our results have here 
typically been obtained for a noise of fixed A = 0.001. 
The number of fingers appearing at early times is related 
to the most unstable wavenumber of the band of unstable 
modes, nevertheless the location and subsequent nonlin- 
ear interaction of the fingers depend on the specific real- 




10 20 30 40 50 60XI0 1 

FIG. 6: Of as a function of time for 15 numerical simula- 
tions realized with the same values of parameters (Pe = 512, 
I = 128, R = 2, s = 1) but different noise r realizations of 
identical amplitude A — 10~ 3 . The dotted and dashed curves 
correspond to the simulations of Figs. and Q respectively. 

ization of the random numbers series. As a consequence, 
it is necessary to compute a set of realizations to get sta- 
tistical information on c/, the main quantity of interest 
here. Figure shows the temporal evolution of aj for 
15 different noise realizations of identical amplitude for 
fixed values of the parameters R, Pe, I and e. As can be 
seen, if fingering starts always at the same onset time t* 
for fixed A, the contribution to the variance due to fin- 
gering saturates to different asymptotic values . This 
corresponds to slightly different nonlinear interactions of 
the fingers as can be seen on Figs. |21 and [7| which show 
the temporal evolution of the fingers for the respectively 
dotted and dashed curves of Fig. If the patterns ob- 
served are very similar during the initial linear phase of 
viscous fingering, the evolution of the fingers is slightly 
different in the nonlinear regime, leading to different val- 
ues of (Too. In particular, merging is observed in Fig. [7| 
leading to the fast development of one finger and, then, 
spreading of the stripe of viscous fluid leading to a larger 
value of doc. 

As a consequence, to understand the influence of fin- 
gering on the broadening of finite slices, it is necessary to 
study the parametric dependence of < Ooo >, the statis- 
tical ensemble averaged asymptotic value of the fingering 
contribution to the variance. 



VI. PARAMETER STUDY 

The quantity < > gives information on the influ- 
ence of viscous fingering on the broadening of finite sam- 
ples. In applications such as chromatography and dis- 
persion of contaminants in aquifers, such a broadening is 
undesirable and it is therefore important to understand 
the optimal values of parameters for which < (7^ > is 
minimum given some constraints. In that respect, let 
us first consider a porous medium in which dispersion is 



FIG. 8: Influence of the sample length I on (Too for (T) Pe = 64 
and (•) Pe = 128 (R = 2, e = 1). 



FIG. 7: Density plots of concentration for the same values of 
parameters and same times as in Fig. |5J but a different noise 
r in the initial condition. 

isotropic (e = 1) and let us analyze the subsequent influ- 
ences of I, R and Pe. The anisotropic case (e ^ 1) will 
eventually be tackled. The mean value < Coo > is plot- 
ted for various values of the parameters, the bar around 
this mean value spanning the range of asymptotic data 
between the minimum and maximum observed. 



A. Influence of the sample length / 

The sample length I has been measured to have practi- 
cally no influence on the onset time t* of the instability. 
Although Nayfeh has shown that the stability of finite 
samples could be affected if the two interfaces are close 
enough [2J], we note that, for the smallest value of sam- 
ple length I considered here (I = 32), the rear interface 
features the same initial pattern as the one appearing on 
the interface between two semi-infinite regions of differ- 
ent viscosities for a same random sequence in the seeding 
noise. Our samples are thus here long enough for the 
onset time t* to depend only on the amplitude A of the 
noise seeding the initial condition and not feel the finite 
extent of the sample. The length I influences nevertheless 
the broadening of the peak and thus < > in partic- 
ular for small /. The points reported in Fig. [S] for two 
different Pe are obtained for one realization and a same 
seeding noise r in the initial condition, leading to a typ- 



ical value of Coo. The smaller the extent I of the sample, 
the sooner the dilution of the more viscous solution into 
the bulk of the eluent and thus the less effective fingering. 
Above a given extent l c , Coo is found to saturate. At first 
sight, this might appear counterintuitive as one could ex- 
pect that, for longer samples, fingering is maintained for 
a longer time thereby enhancing the fingering contribu- 
tion to the variance. A closer inspection to the finger 
dynamics shows on the contrary that, after a transient 
where several fingers appear and interact, only one single 
finger remains (see Figs. 2 and 7). In the absence of tip 
splitting, the stretching of the mixing zone becomes then 
exclusively diffusive as already discussed previously by 
Zimmerman and Homsy |l6l Il7| . This is clearly seen in 
Fig. 5 which shows that the mixing length grows as \ft at 
long times after a linear transient due to fingering. Once 
the asymptotic diffusive regime is reached, the contribu- 
tion of fingering to the broadening of the peak dies out 
and <jf saturates to a^. Above a given critical length l c 
of the sample, the same asymptotic single finger growing 
diffusively is reached before the left and right interfaces 
interact. Hence the same value tXoo is obtained for any 
I > l c - Let us note that the switch from the fingering 
to the diffusive dynamics appears later in time when Pe 
is increased. Indeed larger Pe means more fingers that 
can interact for a longer time before the diffusive regime 
becomes dominant. As a consequence, l c is an increasing 
function of Pe as can be seen on Fig. [H] Further studies 
need to be done to understand the role of e and of pos- 
sible tip splitting occurring for large Pe on the existence 
and value of the critical length l c . 

The fact that the contribution of fingering to the 
broadening of the peak saturates beyond a critical length 
of the sample has important practical consequences for 
chromatography: if fingering is unavoidable, one might 
as well load samples of long extent as the contribution 
of fingering is saturating beyond a given l c . For long 
samples, the efficiency of the process depends then on 
the competition between tr^, and Z 2 / 12, the respective 
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FIG. 9: Onset time t* of the instability for increasing val- FIG. 11: Influence of the Peclet number Pe on < > 
ues of R {Pe = 256, I = 128, e = 1). (— ) Best fit of the (I = 128, R = 2, e = 1). 
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FIG. 10: Influence of the log-mobility ratio R on < (Too > 
for (o) Pe = 128,1 = 128, (•) Pe = 128, Z = 512, (A) Pe = 
256, 1 = 128, (A) Pe = 256, 1 = 512 (e = 1). 



lingering and initial length contributions to the peak's 
variance. We can thus predict that for l c < I < (Joo/\/\2, 
the contribution of fingering is constant and dominates 
the broadening while for / > <j oq /\/Y1, the initial sample 
length becomes the key factor. 



B. Influence of the log-mobility ratio R 

It is easy to foreseen that the larger R, the more im- 
portant the viscous fingering effect |g. First of all, linear 
stability analysis of viscous fingering at the interface be- 
tween two semi-infinite domains |l8| predicts that the 
characteristic growth time of the instability decreases as 
R~ 2 . Although already influenced by the nonlinearitics 
and dependent on the amplitude of the noise, the on- 
set time t* measured in our simulations shows the same 
trend (Fig. |3J). Note that, for very small values of R, 
the onset time becomes very large which explains why, 



for samples of low viscosity, fingering might not be ob- 
served during the transit time across small chromato- 
graphic columns or on small scale contamination zone. 
When R is increased, the viscous fingering contribution 
to peak broadening < > is more important (Fig. I10|) 
with a linear dependence suggesting a power law increase 
for larger R. 



C. Influence of the Peclet number Pe 

The Peclet number Pe is typically experimentally in- 
creased for a given geometry by increasing the injection 
flow rate U. As can be seen on Fig. 1111 < > is found 
to increase linearly with Pe. Fingering induced broad- 
ening can thus be minimized by small carrier velocity U 
as expected. However, the exact influence of the car- 
rier velocity U is difficult to trace because practically, a 
change in U also modifies the dispersion coefficients and 
hence the value of e. In our dimensionless variables, U 
also enters into the characteristic time and length corre- 
sponding respectively to D x /U 2 and D x /U. The concrete 
influence of the carrier velocity is thus more complicated 
to trace in reality. For a fixed injection speed, the Peclet 
number can also be varied by changing the width L y of 
the system. The linear dependence of < > on Pe 
is then related to the fact that in a wider domain, more 
fingers can remain in competition for a longer time so 
that a more active fingering is maintained. This also 
implies that l c is an increasing function of Pe. In chro- 
matographic applications, increasing the diameter of the 
column (i.e. increasing L y here in our model) is thus 
expected to dramatically increase the influence of finger- 
ing in broadening. This explains why fingering really 
becomes an issue for wide contamination zones and in 
preparative chromatography where columns of very large 
diameter (up to one meter) are sometimes constructed. 
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FIG. 12: Onset time t* of the instability for increasing values 
of e for different values of the log-mobility ratio: (V) R — 1, 
(o) R = 2, (A) 7? = 3 (Pe = 128, I = 128). 



FIG. 13: Influence of the ratio between transverse and longi- 
tudinal dispersion e on < cjrx > for (V) R = 1, (o) R = 2, 
(A) R = 3 (Pe = 128, I = 128). Insert: same data on loga- 
rithmic scale for e. 



D. Influence of the ratio of dispersion coefficients e 

Figure IT21 shows the influence of the ratio of dispersion 
coefficients e = D y /D x on the onset time of the instabil- 
ity. As expected from linear stability analysis 01 > de- 
creasing e has a destabilizing effect as fingering appears 
then quicker. This is due to the fact that small trans- 
verse dispersion inhibits the mixing of the solutions and 
favors longitudinal growth of the fingers allowing them to 
survive for a longer time. As a consequence, the less vis- 
cous solution instead of being transversely homogeneous 
invades the more viscous fluid preferably in the longitu- 
dinal direction leading to larger mixing zones and hence 
larger < >. Figure PHfl illustrates that decreasing e 
has a dramatic effect on the broadening of the peak. The 
insert shows the same graphics in logarithmic scale for e. 
< Coo > seems to vary as ln(e) at least for small values 
of e. Peak broadening due to fingering is therefore ex- 
pected to be particularly dramatic for chromatographic 
applications where e ~ 0.04. 

VII. CONCLUSION 

Viscous fingering leads to a mixing between miscible 
fluids of different viscosity. In the case of viscous slices 
of finite extent, fingering is a transient phenomenon be- 
cause the mixing of the two fluids leads to an effective 
decrease of the log-mobility ratio in time. Transient fin- 
gering can nevertheless play an important role because it 
contributes to distortion and broadening of the sample. 
In particular, we have shown that, even if the spreading 
of the sample may look Gaussian at long times because 
dispersion has again become the leading transport phe- 
nomenon, the variance of the peak is larger than expected 
because of fingering at early times. We have demon- 
strated this influence by numerical simulations of viscous 
fingering of miscible finite slices characterizing the onset 



time of the instability t* and the contribution of finger- 
ing to the sample's variance. It is important to note 
that quantitative comparison with experimental data is 
difficult because the exact amplitude of the fingering con- 
tribution to the temporal variation of the peak's variance 
depends both on the amplitude and spatial realization of 
the noise seeding the initial condition which varies from 
one experiment to the other. In this respect, we have 
computed the ensemble averaged asymptotic fingering 
contribution to the peak broadening as a function of the 
four relevant parameters of the problem i.e. the initial 
length of the sample I, the initial log-mobility ratio R, 
the Peclet number Pe and the ratio between transverse 
and longitudinal dispersion coefficients e. The broaden- 
ing of the peaks due to fingering is most important for 
large R and Pe but small e while it saturates above a 
given initial length I of the sample. In chromatographic 
columns for which e ~ 0.04, fingering is thus of crucial 
importance particularly in preparative chromatography 
for which the large diameter of the columns lead to large 
Pe and the high concentration of the samples usually im- 
plies large R. Similarly, for soil contamination, fingering 
will be a major problem in the case of stratified media 
such that e < 1. More work is now needed to explore the 
generalization of this first approach to the case where 
both viscosity and density variations as well as heteroge- 
neous chemistry may interplay as is usually the case in 
the applications analyzed here. 
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